Magnetic-field— induced Luttinger liquid 
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It is shown that a strong magnetic field applied to a bulk metal induces a Luttinger- liquid phase. 
This phase is characterized by the zero-bias anomaly in tunneling: the tunneling conductance scales 
as a power-law of voltage or temperature. The tunneling exponent increases with the magnetic 
field as BlxiB. The zero-bias anomaly is most pronounced for tunneling with the field applied 
perpendicular to the plane of the tunneling junction. 

PACS numbers: 71.10Pm,72.15Gd,72.15 
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A strong magnetic field applied to a bulk metal tends 
to reduce the effective dimensionality of charge carriers 
from 3D to ID. This feature is most pronounced in the 
ultra-quantum limit, when only the lowest Landau level 
remains populated. The reduction of the effective dimen- 
sionality in a system of interacting particles is expected 
to result in a number of unusual phases, which are one- 
dimensional irunature, such as spin-and charge-density 
waves (CDW)EjH, Wigner crystaJa, excitonic insulator 
(EI)El, re-entrant superconductor^, etc. A field-induced 
CDW is believed to have been observed in graphite (see, 
e.g., Ref. [ [?J for an extensive bibliography and discus- 
sion). There have been also earlier reports of tentative 
field-induced EI transitions in Bii-^SbJj and BjlL 

In this paper, we focus on another field-induced state, 
which is not related to any instability, but evolves adi- 
abatically from the conventional three-dimensional (3D) 
Fermi-liquid as the magnetic field increases. This state 
is a Luttinger Liquid (LL), whose existence can be antic- 
ipated from the following simplified picture. In a strong 
magnetic field, electron trajectories are helices spiraling 
around the field lines. A bundle of such trajectories with 
a common center of orbit can be viewed as a ID conduc- 
tor ( "wire" ) with the Fermi velocity vp 1 determined from 
the condition of the fixed carrier density n in the bulk, 



Vp = 2n 2 nl 2 B /m\\. 



(1) 



Here Ib — (eB) ' is the magnetic length, and Tom is 
the effective mass along the field; we assume also full spin 
polarization of electrons, and use the units with K = c= 1. 
In the presence of electron-electron interactions, each 
"wire" , considered separately, is in the LL-state. Inter- 
actions with small-momentum transfers among electrons 
on different "wiees" do not change the LL-nature of a 
single- wire stataHJ. 

The quantity of primary interest of this paper is the ex- 
perimentally measurable density of states (DOS) at the 
sample boundary. We calculate this quantity first by 
treating perturbatively the interaction between electrons 
moving in 3D space in the presence of a quantizing mag- 
netic field. We demonstrate that the boundary DOS ex- 



hibits a characteristic for a Luttinger liquicfl power-law 
anomaly at the Fermi level: 



v(s) cx 



e — e. 
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where ep 1 = mil (vp) 2 /2 is the Fermi energy of the ID 
motion along the magnetic field. 

We establish the correspondence between a 3D electron 
plasma in a quantizing magnetic field and a Luttinger 
liquid by considering the electron-electron interaction in 
the basis of coherent states (CS) on the von Neumann 
latticed 2 ]. In this basis, a system of 3D electrons in the 
ultra-quantum limit is equivalent to a lattice of ID wires 
pointing in the direction of the field and separated by 
distance v / 2rfs, similar to the simplified picture men- 
tioned aboveEJ. The problem of an array of wires coupled 
via the long-range Coulomb interaction [Dzyaloshinskii- 
Larkin (DL) modelEH] can be. solved exactly .either via 
summation of diagram serieaij or bosonizatiorEj. Utiliz- 
ing the latter technique, we derive the "bulk" tunneling 
density of states. This quantity also exhibits a power-law 
anomaly (Q) with a different exponent a. Different values 
of bulk and boundary exponents is another characteris- 
tic feature of a LLtl. We calculate explicitly the values 
of the exponents a and a s in the limit of a long-range 
interaction potential, 



i B e B < i, 



(3) 



where k b = uip/vp is the (field-dependent) inverse 
screening length, and lo v is the zero-field plasma fre- 
quency. Condition (||) is satisfied in magnetic fields 
B < Bq = ir en/m\\e 3 , where e is the background dielec- 
tric constant. The introduced here characteristic mag- 
netic field Bo can be expressed as Bq ~ B q /r s in terms 
of the magnetic field B q de-populating all but the low- 
est Landau level, and of the gas parameter r s . Under 
condition ([s|), we find 
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The tunneling current-voltage characteristic is thus a 
power-law 



1 



I oc |y| as+1 sgnV (5) 

for e\V\ <C Ep. Eq. (||) refers to the situation when 
the magnetic field is perpendicular to the plane of the 
tunneling contact. If the field is parallel to the contact 
plane, electrons move along skipping orbits. Interactions 
between electrons moving in the same direction do not 
result in anomalous scaling, therefore there is no tunnel- 
ing anomaly for this field orientation. This means that 
by tilting the field one can vary the tunneling exponent 
from its maximum value, given by Eq. (Q), to zero. 

Another indication of the tunneling anomaly comes 
from the gapless behavior of a plasmon in a strong mag- 
netic field. For uj c 3> u>, where lu c = eB/m c is the cy- 
clotron frequency, the dispersion relation for a classical 
magnetoplasmon is given by u 2 = (k 2 /k 2 }u 2 [ |l6). As we 
see, >0 for k z ^0 at finite k. A gapless charge mode 
slows down the relaxation of an excessive charge added 
to the conductor in a tunneling event. This mechanism is 
responsible for the zero-bias anomalies in tunneling into 
disordered metals and ID Luttinger liquids. A metal in a 
strong magnetic field provides one more example of such 
a behavior. 

The LL is destroyed either by backscattering between 
electrons on different "wires", which results in a CDW- 
instability, or by formation of bound electron-hole pairs 
(for the case of a semimetal), which leads to an EI- 
instability. Also, impurity scattering transfers electrons 
between the "wires" and thus destroys the ID motion. 
The relevant energy scales for these processes is the en- 
ergy gap, A, of either CDW- or El-origin, and the level 
broadening, T, respectively. On the other hand, a LL- 
behavior sets in at energies smaller than the Fermi energy 
Ep . The gap A increases with the field, whereas Ep de- 
creases with the field. The LL-state should thus exist in 
the energy interval which narrows down as B increases: 

max{A(B),r} < e < ef oc B~ 2 . (6) 

Before giving the derivation of the results announced 
above, we would like to discuss their relevance to the ex- 
periment. The search for LL-like tunneling anomalies in 
truly ID. systems (edges of 2D electron gas in the jEDHE 
rcgimeO, carbon nanotubesES, and quantum wiresES) has 
been a very active but remarkably difficult area over the 
last few years. Compared to experiments on tunneling 
into ID systems, tunneling into a 3D system in the ultra- 
quantum limit has the obvious advantage of the macro- 
scopic system size. It is also advantageous that the tun- 
neling exponent is a function of the external parameters 
(magnitude and direction of the magnetic field), which 
can be varied over a wide range. 

The choice of the right material is crucial. Among the 
low-density semimetals (Bi and its alloys and graphite) 
used conventionally for high-field studies, graphite seems 
to be the optimal candidate. For B along the c-axis, 
the ultra-quantum limit is achiej^ejd at B ~ 7 T. With 
material parameters for graphite£3 (m|| — IOtoq, e = 6, 



and n = 3x 10 18 cm~ 3 ), Eq. (Q) gives a = 0.8 already at 
7T. For comparison, an anomalously large value, of the 
background dielectric constant (e ~ 100) for BE!l and a 
smaller value of m|| (= 0.5mo for holes) makes a to be 
very small: even for i? = 100.T, one has a = 0.05. Accord- 
ing to a recent experiments 2 ], the field-induced CDW gap 
in graphite A < 1 K for B < 30 T and takes its maximum 
value (A ~ 10 K) at B ~ 50 T. For fields in the range 
7T <£> <30T, inequality (jfj) is satisfied for more than 
two decades of energies (voltages) . 

In the limit of weak electron-electron interactions, one 
can derive the field-induced zero-bias tunneling anomaly 
(H) by calculating interaction corrections to the trans- 
mission coefficient through a tunneling barriered. These 
corrections can be then summed up by using the renor- 
malization group procedure. 

Let a potential barrier of transmission coefficient to <C 1 
separate two metallic half-spaces, z < and z > 0. The 
magnetic field is perpendicular to the contact plane. We 
assume that the base-electrode (z < 0) is made of a high- 
density metal and therefore is not affected by the mag- 
netic field. Electron-electron interactions are then not 
important in this half-space and we treat it as a Fermi 
gas. The counter-electrode (z > 0) is a low-density metal 
in the ultra-quantum magnetic limit. In this half-space, 
electrons interact via potential U(r — r'). Keeping only 
the leading terms in to<^ 1, the free wave-function in the 
Landau gauge A = (— yB, 0,0) for z > and far away 
from the barrier is given by 

K*,-p* = ~ ( 2l siiw/yZT) (r ± ) (7a) 

C^^'V^tW, (7b) 

where p z = \j2me > 0, rj^ = (x, y) and 

<(r±) = (V^sL.r^V^exp [-(y + pj%) 2 /2£ 2 B ] . 

The leading correction to the transmission coefficient 
comes from the exchange interaction^, if the interaction 
potential is smooth on the scale of Fermi wavelength, see 
Eq. (|). 

The correction to the wave caused by the exchange 
potential V x has the form 

*W«=*L P >) + / J dv'dv"Gl(r,r') 

x^(r',r")$° t!fc (r"), (8) 

where Gf (r,r') is the Green's function of the free 
Schrodinger equation in the half-space z > 0, and e z = 
p z /2m 2 is the energy of electron motion along the field. 
The integration in (|^) goes only over those regions where 
electrons interact, i.e., z',z">0. 

The renormalized transmission coefficient is extracted 
from the asymptotic form of f PijPz for z — > +oo which, 
in its turn, is determined by the asymptotic behav- 
ior of G^(r,r') [see Eq. (||)]. The asymptotic form of 
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the Green's function in the infinite space can be writ- 
ten as G Pi (r,r')| z ^ +00 = G||(z,z',p z )G_ L (r J _,r^ L ), where 

G'e z (z, z') — (m\\/ip z )e lPz \ z ~ z ' is the Green's function for 
the ID motion along the z-axis, and 



G±(r ± ,r' ± ) 



2tt£ 2 b 



-[r ± -r' ± ] 2 /4e 2 B -i(y+y')(x-x')/2e 2 E 



(9) 



The asymptotic form of Gf(r,r') is then obtained via 
the method of images 



The exchange potential is given by 



V x (v,v')=-U(v-t 



') f Fd l± f° r*o fr01 * 

>J 2n 2n 1 >l 



_ U(T _ ^f*+/) G±{tXtt , ±)i ( 10 ) 
tt[z + z') 



where p B = m\\Vp. In Eq. (|1(J), we retained only that 
term which at p z « p B leads to a logarithmic divergence 
in the integral over z + z' in Eq. (||). 

The transmission coefficient in the first order with re- 
spect to interactions is 



t = tn 1 



In ■ 
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k b Ib~ i B \p,-pf 
where the inverse screening radius 

k b = Wpfvf = 2od B ~ 1 oc B. 



In- 



1 



(11) 



(12) 



This result is valid if both logarithmic factors are large. 
The contribution t from the Hartree term is smaller than 
( pi] ) by \a(\/ k b £b)- The higher order terms in the ex- 
pansion of t can-be summed up via the renormalization 
group procedures. The result is 



t oc \p z 



p3\ 



,/2 



oc |e 
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(13) 



As I oc \t\ 2 V, Eq. (|l3l) yields Eq. (§) with exponent (|). 
The perturbative result (|l^) is valid for a <C 1, which 
means that a s is also small within this approach. 

We now demonstrate how the LL-state emerges in the 
coherent state formulation. The method of bosonization, 
applied in this formulation, provides a siaaple tool for 
calculation the "bulk" ID density of statesJlil. Also, it al- 
lows us to estimate the tunneling exponents for stronger 
interaction (a > 1). The ^-operator is expanded over 
the coherent states basis 



(14) 



N R 



where r=(z, rj_), H=(R X , Ry). A coherent state xjvr.(i"_l) 
is a simultaneous eigenstate of the energy operator 



(2m.)- 1 (p ± - eA) 2 X nr = u c {N + 1/2) X nr, 
and of the "guiding center" operator 

(x - i yo)XNTL = {Rx - iRy)XNR, 

where xq = x + (muj c )- 1 (py — eA y ) and j/q = y — 
(mui c )~ 1 (p x — e.A x ) are the operators corresponding to 
the classical coordinates of guiding centersa. In the sym- 
metric gauge, A = (1/2)B x_c. an explicit form of the 
CS corresponding to N = 



G>(r,r')U+oo = \GUz,z')-Gl(z,-z') \ G x (t x ,t' x ). XoR (r ± ) 



2tt£, 



■ exp 



(r ± - R) + 2ir ± A R 
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Cl x by 



dy^x ■ 



where a A b 

As any set of coherent states, the set of xnr remains 
(over) complete even if defined on a subset of points in the 
complex plane, R x + iR y , rather than on the whole plane. 
The von Neumann theoremEa allows one to choose this 
subset as sites of a 2D square lattice. To fix the value of 
the lattice spacing, one notices that a total of n particles 
per unit volume are now distributed over "wires" , which 
are parallel to the magnetic field and cross the plane per- 
pendicular to the field at the sites of the von Neumann 
lattice. For N = 0, the number density per unit length 
of each "wire" is Pp/n, their the areal density is 1/s 2 , 
thus n — p B /its 2 . Comparing this equation to (0), one 
gets s — \/2tt£ b ■ This is a maximum value of the lattice 
spacing for which the set is still complete. States of the 
CS basis on the von Neumann lattice describe particles 
which are localized on the lattice sites in the direction 
transverse to the field but free to move along the field. 
This picture is in a close resemblance to the classical tra- 
jectories spiraling around the field lines. 

Coherent states are not orthogonal. In particular, 



(XorIxor') = g ex P 



(R R) 



2iR A R 



U 2 B 



(15) 



which can be replaced by 2n£ 2 B 5 2 {YL — R') in the limit 
£b — * 0. Only in this limit the coherent set diagonalizcs 
the quadratic part of the Hamiltonian. Physically, this 
limit means that £b is much less than the characteris- 
tic value of |R — R'|. For interacting electrons, typical 
|R — R'| is of the order of the screening radius 1/k b . 
Thus the approach based on coherent states is valid for 
a long-range interaction as defined by condition (Q). 

Using asymptotic orthogonality of the coherent states, 
one arrives at the effective ID Hamiltonian (from now on 
we concentrate on N = and suppress index N) 



H = 



R 



dzaj^(z) 



2 



Pz 



2mi 



or (2) 



- y 

2e ^ 

R.R' 



dzdz 



,a f R (z)al l ,(z')a RI/ (z')a R (z) 



(R - R') 2 + (z - z') 2 



(16) 
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Due to the long-range nature of the Coulomb potential, 
it suffices to keep only the forward scattering processes in 
(H) (both for R = R' ancLR ^ R), in which case (|l6|) is 
identical to the DL modelELpBackscattering for R^R' 
leads to a CDW instabihtyeliI3. The forward scattering 
part of ( |l6|) is diagonalized via bosonizing the fermions 



-i( P fz+^{<f>R+e R }) + e i( P fz+^{<t> R -e R }) . 



where [4> R (zi), d Z2 9 RI (z 2 )} = «<5r,r/<5(zi - z 2 ). The re- 
sult of a straightforward calculation for the equal-point 
(Matsubara) Green's function is 

0(r) ~ T - l e -K T \ 



f(r) = 2^4 



2 > I d2f l±- [ Ad( lz ( u - l) 2 1 - e~ uv Fi* T 



BZ 



(2tt) 2 J q 2tt 



u 2 ^l + (k b ) 2 £(Q + G)- 2 , Q = (qx, q,). 

G 

Here the ^-integration goes over the Brillouin zone (BZ) 
of the von Neumann lattice, G are the reciprocal lattice 
vectors, A is the cutoff, and uv B has the meaning of a 
plasmon velocity. Typical q z in /(r) are determined by 
the low-energy scale of the problem (max{T,eV}) and 
thus small, which allows one to neglect q z in function u. 
The DOS reduces to the form (0) with an exponent 



a 



nii 



BZ 



d 2 qi_ (u - l) 2 
(2tt) 2 u ' 



(18) 



For k b £b^ 1, the sum over the reciprocal vectors in (118) 
is dominated by the G = term so that u 2 ~ l+(n B /q±T ■ 
The main contribution to the integral over q± in (|18]) 



comes from the region q± 



C 1/ffi - G, thus the 



integration can be extended over the entire g-space. The 
bulk DOS behaves as |e— e B \ a , where a is given in Eqjfl). 

Note that in the main contribution to a [see Eq. (pj|)], 
a typical value of u — 1 is of the order of unity. This 
indicates a non-perturbative nature of the result for a. 
The lowest-order perturbation theory, applied to the bulk 
case, would have given a cx e 4 , whereas the correct result 
is a tx e 2 . 

Tunneling into the end of a semi-infinite sample can 
be treated by imposing the boundary condition on the 
current operator: j R = (l/i0r)<9 r </>R = (1/ s/tt)8 x 6 r = 
for z = 0. This translates into the boundary conditions 
for the propagators: V r (z,z',t) — (p R {z, r)p R (z', 0) — 
p R {z, 0)), where p = <f),9. V R can be constructed from 
the propagators for the infinite medium, P R , by using 

the method of images. Evidently, V R — at the bound- 
ary whereas V R (0, 0, r) = vr" 1 / dq z P R (q z ,T). The DOS 
at the boundary is again of form (^|) but with an exponent 



d 2 q± 



BZ 
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(19) 



In contrast to the bulk case, the logarithmically divergent 
integral over q± in (pit) is cut off at the BZ boundary, i.e., 



for q± ~ 1/Ib- Therefore, the condition of asymptotic 
orthogonality of CS is satisfied only with the logarithmic 
accuracy. The result for a s is given by Eq. (|]). 

In the strong-coupling case (a > 1), typical q± ~ 1/^b, 
which means that the approach based on CS becomes 
invalid quantitatively. However, one can still estimate 
tunneling exponents (both for the "bulk" and "edge tun- 



neling) as a — Cye 2 /vp, where the numerical constant 
C ~ 1 cannot be determined by this method. 
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